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Abstract 

We present theoretical and experimental data on the threshold behaviour of nanodot creation 
with swift heavy ions. A model calculation based on a two-temperature model taking the spatially 
resolved electron density into account gives a threshold of 12 keV/nm below which the energy 
density at the end of the track is no longer high enough to melt the material. In the correspond- 
ing experiments we irradiated SrTiC>3 surfaces under grazing incidence with swift heavy ions. The 
resulting chains of nanodots were analyzed by atomic force microscopy. In addition, samples irradi- 
ated under normal incidence were analyzed by transmission electron microscopy. Both experiments 
show two thresholds, connected to the appearance of tracks and to the creation of fully developed 
tracks, respectively. The threshold values are similar for surface and bulk tracks, suggesting that 
the same processes occur at glancing and normal incidence. The experimental threshold for the 
formation of fully developed tracks compares well to the value obtained by the theoretical descrip- 
tion. 

PACS numbers: 68.37.Ps, 68.37.Lp, 61.80.Jh, 61.82.Ms 
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I. INTRODUCTION 



The irradiation of solid matter with heavy ions of MeV energy has long been known to 
result in structural modifications ranging from stoichiometric or geometric defects and amor- 
phization in the bulk up to the creation of hillocks and chains of hillocks on the surface^ 1 2 1 3 1 4 1 5 1 . 
In this energy range the modification of the material is not due to direct collisions of the 
projectile ion with the atoms of the target material (nuclear stopping power regime) but 
rather to a very intense interaction of the projectile with the electronic system of the tar- 
get (electronic stopping power regime). How this electronic excitation is transformed into 
material modifications depends very much on the physical properties of the target material 
itself. Most insulating materials exhibit nanosized hillocks at the impact zone of the ion 
when irradiated normal to the surface. This effect has been known and studied in detail for 
quite some time (see e.g."™) but is still not fully understood. 

A rather successfull approach to describe the transformation of an electronic excitation 
into a heated lattice is based on a two-temperature model (TTM) and requires solving the 
coupled differential equations for the electron and the phonon system, respectively. The 
model has been used to explain e.g. track radii in various materials created by irradiation 
perpendicular with respect to the surface under the condition that the electron-phonon 
coupling constant g is a fitting parameter. 

Recently, the two-temperature model has been modified to include the spatially resolved 
electron density instead of a homogeneous free electron gas^ which is important if the irra- 
diation takes place under glancing angles. In this case, a single ion is able to create a chain 
of hillocks along its otherwise latent track within the bulk. According to the model, the 
chains occur due to the non-homogeneous nature of the electron density. Every time the 
projectile travels through a region with a high density the energy loss is sufficiently high to 
feed energy into the electronic system. This energy is transferred into the phonon system, 
a local melting occurs which finally results in a nanodot on the surface. Since the electron 
density corresponds to the periodicity of the crystal, the nanodots appear with a certain 
periodicity on the surface. 

So far, these chains of hillocks have always been produced with ions experiencing a 
stopping power of about 20 keV/nm and nothing is known about the morphology at lower 
stopping powers. The aim of this paper is to study ion induced nanodot chains in SrTiC>3 as 
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a function of kinetic energy of the projectile in order to determine a threshold value for the 
stopping power. In addition we wish to show that the recently presented theoretical approach 
can be applied to predict a threshold value which is in agreement with experimental data. 

This paper is organized as follows: At first we will present experimental data from two 
different methods to determine the threshold energy for the production of nanodots: In 
section A we discuss the use of atomic force microscopy (AFM) to determine the minimum 
stopping power required to create a chain of nanodots on the surface. As we normally treat 
the surface chains as tracks on the surface, it is interesting to compare the two manifestations 
of a swift-ion passage: surface chains and tracks in the bulk. Therefore, we present in section 
B results obtained with transmission electron microscopy (TEM) in order to determine the 
stopping power required for amorphization. We then derive a value for the threshold from 
theory and discuss the experimental data with respect to our theoretical predictions. 

II. EXPERIMENTAL DETERMINATION OF THRESHOLD VALUES 
A. Grazing incidence 

The single crystal samples of SrTiO3(100) (Crystec, Berlin and Kelpin, Neuhausen) have 
been irradiated without prior surface treatment at the 6.0 MV Tandem Van de Graaff ac- 
celerator at the Institute Ruder Boskovic in Croatia^. The irradiation was done using I 
ions with different kinetic energies resulting in different stopping powers (calculated with 
SRIMP): 6.5, 13, 18, 23, and 28 MeV yielding 2.9, 5.3, 7.2, 9.0, and 10.5 keV/nm, respec- 
tively. The angle of incidence with respect to the surface was kept at <ft = 1.3°. The samples 
were oriented along the (001) direction within a few degrees. Fluences were typically chosen 
to yield around 10 tracks per /zm 2 on average, assuming that every ion produces one track. 
For lower stopping powers also higher fluences were used to compensate for a possible loss 
in production efficiency Additional samples of SrTiO^lOO) were irradiated under = 1° 
at the beam lines IRRSUD (92 MeV Xe 23+ , 21 keV/nm) and SME (700 MeV Xe 23+ , 29 
keV/nm) at the GANIL in Caen, France. After irradiation the samples were analyzed by 
AFM in the tapping mode under ambient conditions. Some samples have been cleaned prior 
to AFM-measurements by snow-cleaning and wiping with ethanol. By comparing chains of 
nanodots from irradiated samples with and without cleaning treatment we made sure that 
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Figure 1: Topography images of a SrTiOs surface irradiated with different energy losses under 
grazing incidence. I ions under 1.3° were used in a)-c), Xe ions under 1° in d) and e). Frame size: 
1 /im x 1 //m. To enhance the contrast false colouring was used. Electronic energy losses in the 
sub-figures are a) 5.3 keV/nm, b) 7.2 keV/nm, c) 10.5 keV/nm with a height profile of the marked 
track, d) 21 keV/nm and e) 29 keV/nm, respectively. 

this procedure does not affect neither the shape nor height of the hillocks. 

Fig. shows an AFM image of the typical chains that are produced on a SrTiC>3 surface 
after irradiation with 10.5 keV/nm ions at a glancing angle of <p — 1-3°. Each chain in fig. ^ 
is produced by a single ion because the average number of chains (~ 10 per /im 2 ) corresponds 
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well with the nominal fluence of lxlO 9 ions/cm 2 . Thus, the production efficiency is close 
to one. The chain length is about 300 nm and the individual hillocks within the chains are 
a few nanometers high. At this stopping power height and diameter of the dots appear to 
be the same as found in earlier experiments with higher stopping powers, see figs. [l]i and 
For a general discussion of the velocity effect in the electronic energy loss region on the 
measured chains of hillocks, sedPl 

In fig. [pi andfij) we show a representative AFM image of samples that have been irradiated 
with ions of 5.3 keV/nm and 7.2 keV/nm, respectively (same fluence, same angle as in fig.fTj^). 
Ions with 5.3 keV/nm (fig. [T£i) as well as with lower energy losses produce no visible chains. 
First chains were detected at 7.2 keV/nm (fig. fib). As can be seen, the average chain length 
as well as the number of chains is smaller than for higher energy losses. Our statistic was 
not sufficient to quantitatively discuss the chain production efficiency (chains per incident 
ion), but we were able to study the length of the chains. 

To demonstrate the evolution of chains as a function of ion energy more clearly we have 
performed a statistical analysis of at least 60 chains for each kinetic energy, see fig. [2] It 
should be noted that the error bar in the length measurement does not originate from uncer- 
tainties in the measurement itself, but from the statistical nature of the interaction process. 
From fig. [2]we find that the minimum stopping power to create a chain is between 5.5 keV/nm 
and 7.5 keV/nm. Up to a stopping power of between 10.5 keV/nm and 21 keV/nm the chain 
length increases with increasing kinetic energy until a constant chain length is reached. The 
SRIM calculations as well as our calculations in section III were made for Xe ions. Because 
iodine is right next to xenon in the periodic table, there is not much difference between those 
two ions with respect to the modifications they produce. Experimentally, we have seen no 
difference in chains created by iodine and Xe ions. 

B. Perpendicular incidence 

A SrTiC>3 single crystal (Crystec Berlin), cut along the (lOO)-plane, was irradiated under 
normal incidence with 92 MeV Xe 23+ ions with a fluence of 2x 10 13 ions/cm 2 at the IRRSUD 
beamline at GANIL, Caen, France. At these fluences, individual tracks are covering each 
other. After irradiation the sample was cut and glued together with epoxy, in order to 
produce two layers facing each other. This sandwich was then encased and glued in a 3 mm 
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Figure 2: Chain length as a function of stopping power. The data points at 21 keV/nm and at 
29 keV/nm are taken from samples irradiated with Xe ions under 1°, the chain length was corrected 
to present an incidence angle of 1.3°, i.e. they were multiplied with cot 1.3°/ cot 1° = 1.3. The other 
data points are from irradiations with iodine ions under 1.3°. 

thin walled tube. The tube was cut into 1 mm thick discs. These discs were mechanically 
thinned to about 100 /im in thickness and dimpled so that the center thickness was less than 
5 jum. Finally, a Gatan's Precision Ion Milling system was used to thin the specimen from 
both sides with 3.5 keV Ar + ions at 5° incidence angle until perforation. These targets were 
studied with a 2010F Jeol Transmission Electron Microscope equipped with a 200 keV field 
electron gun. 

A TEM-image of that cross section is shown in the upper part of Fig. |3j The bright 
triangle on the right side stems from the glue, denoting the surface region of the crystal. 
The ions penetrate the crystal in the direction indicated by the arrow. In the first part 
along the trajectory, the crystal was amorphized by the projectile, as can be seen from the 
diffraction pattern on the right side, taken from this part of the crystal. This first part 
extends from the surface to a penetration depth of ~ 3.6 /im. To correlate this data with 
the stopping power we performed SRIM calculations. The result is shown in fig. |4} In 
the first part, where the crystal was amorphized, the electronic energy loss decreases from 
21 keV/nm to 11.7 keV/nm. From this we determine the threshold of amorphization to be 
about 11.7 keV/nm. 

At larger penetration depths, partial amorphization is seen in fig. [3] This continues to a 
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Figure 3: Top part: TEM image of a SrTi03 sample irradiated under normal incidence. Ions enter 
on the right side, see arrow. The axis below the image shows the depth from the surface. Lower 
part: center image is a magnification of the central part of the image on top. The left and right 
images on the lower part show the diffraction patterns obtained at ion penetration depths of 6 /^m 
and 3 fim, respectively. 

depth of 5.1 /mi. At this depth, the energy loss is 7.3 keV/nm, according to SRIM, shown 
in fig. |4j At even larger penetration depths, the material is still mono-crystalline, as seen by 
the diffraction pattern on the left in fig. |3j Therefore we can distinguish two thresholds, one 
at 7.3 keV/nm for the appearance of partially amorphized tracks and one at 11.7 keV/nm 
for complete amorphization. 

III. THE TWO-TEMPERATURE MODEL 

After we have experimentally determined the minimum stopping power required to create 
a chain of nanodots we focus on the theoretical description. Unfortunately, theoretical 
models describing the hillock formation in detail are still lacking. However, because we are 
only interested in the threshold energy a modified two-temperature model can be used. The 
basic procedure to determine the threshold energy is described in detail inP', and will be 
explained here only insofar, as is necessary to understand the additions we have made to 
the model. 

In order to determine the minimal necessary stopping power for hillock creation the 
following scheme was applied. The calculations were performed on the basis of a 3D two- 
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Figure 4: SRIM simulation of Xe ions impinging on a SrTi03 single-crystal at 92 MeV kinetic 
energy. Shown is the electronic energy loss as a function of the penetration depth of the incident 
ion. Thin lines denote the thresholds for complete amorphization and for the appearance of partially 
amorphized tracks, respectively, as determined from the TEM images in fig. 3. 

temperature model using space- and time-resolved electronic stopping powers for the actual 
case. It was proven before, that there is a fixed maximum depth for the incident ion from 
where it is possible to create surface hillocks due to the electronic energy los ^ 12 l Therefore, 
the travelled length of the incident ion at the point perpendicular to the last hillock can 
easily be extracted from the incident angle and the above-mentioned maximum depth out of 
geometrical calculations. In the next step, we calculate the energy loss of the ion at exactly 
this point for different kinetic energies and use it as a source for the 3D TTM. In this way, 
surface temperatures at the end of the hillock chain can be obtained as a function of the 
incident ion energy, see fig. [5} Assuming that the melting temperature {Tmeit = 2353 KpJ) 
has to be reached in order to induce surface modifications, the incident ion has to have a 
minimum kinetic energy of 31 MeV, corresponding to an energy loss of 12.7 keV/nm. This 
calculation ignores the necessary (but unknown) heat of fusion, and the energy loss obtained 
in this way therefore presents a lower limit. Additionally, this model calculation was done 
for Xe ions, the minimum kinetic energy value for I ions would be slightly higher to achieve 
the same stopping power. 
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Figure 5: Temperature in K at the surface of SrTi03 if irradiated by Xe ions under different 
kinetic energies. To reach the melting temperature of SrTiOs the projectile has to have an energy 
of 31 MeV corresponding to a stopping power of dE/dx =12 keV/nm. The thin line represents a 
linear fit to the data. 



Grazing angle 


Appearance of hillocks Constant length of chains 




5.3 - 7.2 


10.5 - 21 


Normal incidence 


Appearance of tracks 


Amorphization 




7.3 


11.7 


Theory (grazing angle) 




12.7 



Table I: Threshold values (in keV/nm) obtained by experiments and theory. 



IV. DISCUSSION 

In table I, all the threshold values obtained for glancing and normal incidence, as well as 
the theoretical prediction are shown. 

Under glancing angles, tracks show up at kinetic energy losses larger than (5.3 - 
7.2) keV/nm, when studied by means of AFM (section IIA). This fits well with the TEM 
measurements (section IIB), where the appearance of tracks was seen for energy losses of 
7.3 keV/nm. 

The TEM measurements show at lower energy loss values discontinuous tracks; and the 
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AFM measurements reveal for lower energy losses only tracks which are shorter than those 
with higher energy losses. Amorphization in the bulk as well as fully developed surface tracks 
can only be seen at energy losses above 11.7 keV/nm and (10.5 - 21) keV/nm, respectively. 

This second threshold fits well with the calculated threshold of 12.7 keV/nm, where 
melting of the surface over the full surface track length was defined as a requirement in 
order to produce hillocks, see section III. As was shown before^ the track length for high 
energy losses is given by a maximum depth d, from which surface modifications are possible. 
The depth is for SrTiC>3 d ~ 9 nm, and the track length / can then be calculated by 
/ = d cot <p. For sufficiently high energy losses, d is constant even if the initial energy loss of 
the impinging ion changes by over 709$^. If the energy loss approaches the threshold, the 
remaining energy loss at the maximum depth d may not be enough for surface modification. 
This can explain, why we observe increasing surface track lengths with increasing energy 
loss in fig. 2. 

The theoretically predicted threshold assumes that the total chain length is constant in 
order to make assumptions on the necessary stopping power. The calculation of the lower 
threshold is in principle possible as well, but this would require an additional variation of 
the surface track length which would tremendously increase the computational effort. The 
discovery of two different values for the hillock creation is not a contradiction, it merely 
reflects that the hillock creation at stopping powers > 12.7 keV/nm are independant of 
statistical fluctuations. Below this value the individual hillock creation is a more statistical 
phenomenon which results in "incomplete" tracks. 

The determination of the threshold with chains has two main advantages. Experiments 
performed under perpendicular incidence create single hillocks which are difficult to identify 
in ex situ experiments, especially if the production efficiency goes down or the single hillocks 
become smaller when approaching the threshold. A chain is easily spotted and can be 
unambiguously identified as an ion-induced feature, for example the direction must comply 
with the direction of the irradiation. Since we made no specific assumptions with respect 
to the material, we believe that the method is capable of predicting the minimal energy 
required to create nanodots in many other insulating materials of technological importance, 
such as AI2O3 or TiC>2. The determination of the threshold energy by means of analyzing 
surface tracks has proven useful even for systems like SiOj-^or III-V semiconductors^ where 
the connection to the electron density might not be as straightforward. 
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